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Abstract 



I review recent determinations of the (on-shell) charm- and beauty-quark masses in 
the framework of relativistic and non-relativistic ratios of Laplace transform QCD 
^ ' moment sum rules. The validity of the non-relativistic version of QCD sum rules in 

' this particular application is discussed. 
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With the advent of the Heavy Quark Effective Theory (HQET) there has been a revived 
interest in apphcations of QCD sum rules to the heavy quark sector. The purpose being 
the update of old determinations, as well as the performance of new calculations to ex- 
tract accurate values of various dynamical quantities entering the HQET. Chief of these 
quantities is the ubiquitous heavy-quark mass. I review here recent determinations of 
the (on-shell) charm- and beauty-quark masses performed by confronting very accu- 
rate experimental data on the charmonium and the upsilon systems with ratios of 
relativistic and non-relativistic Laplace transform QCD moments. The latter theoretical 
framework, suggested by Bertlmann 0, offers several advantages, e.g. radiative and non- 
perturbative corrections are well under control, and the non-relativistic limit follows quite 
naturally from quantum mechanical analogues @]. This version of QCD sum rules leads 
to an expansion in powers of the inverse of the heavy quark mass which allows one to test 
the range of validity of the non-relativistic limit, and more generally, to assess the role 
of mass corrections. This might be of interest for calculations based on the simplifying 
assumption ^qcd/^q ^ 1- Also, a comparison of the results from the relativistic and 
the non-relativistic determinations provides an estimate of the systematic uncertainties 
affecting this technique. 



I begin by considering the two-point function 

n^.(g) = t J d^xexp{tqx){0\T{V^{x)V,+ m\0) = {-g^,q^ + g^g.)H(g2) , (1) 

with Vf^{x) = Q{x)'~^^Q{x). The function H(g^) has been calculated in perturbative QCD 
at the two-loop level H], with its imaginary part given by 
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where f = y 1 — 4mQ/s, and mq is the charm- or the beauty-quark on-shell mass: mq = 
"^q(Q^ = "^q)- The leading non-perturbative term in the operator product expansion of 
H(g^) involves the gluon condensate, i.e. 
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The function H(g^) satisfies a once-subtracted dispersion relation, and the subtraction 
constant can be disposed of e.g. by taking the Laplace transform 
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The quantity of interest to us here is the ratio of the first two Laplace moments, which 
can be expressed as 

7^(a) = — ^ In H(a) . (5) 
da 
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Substituting Eq.(2) into Eq.(4) one can carry out the integration analytically. The result 
for the ratio Eq.(5), with ui = Amqa, is 
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The function G{b,c,uj) is related to the Whittaker function Wx^fj,{uj) through p 

G{b,c,u)=u;>'-'/^e^/^Wx,,{u;) , 
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with n = {c- 6)/2, and A = (1 - c - h)/2. 



The above expressions involve no approximations, other than the two-loop perturbative 
expansion, and the truncation of the operator product expansion beyond the leading non- 
perturbative term. We shall refer to Eq.(6) as the fully relativistic Laplace ratio. In the 
non-relativistic (heavy quark-mass) limit, the Laplace transform Eq.(4) becomes 

roo 

n(r) = / dEe-x.Y>i-TE)lmTl(E) , (18) 
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where r — Amga, and s — {2mQ + E)"^ so that £^ > 0. The Laplace ratio Eq.(5) is now 
given by 

7^(T) = 2mQ - In n(r) . (19) 
or 

After expanding the functions G{b,c,u!) entering Eq.(6), one obtains the non-relativistic 
ratio 
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where the expansion has been truncated at the next-to-next to leading order in I/itlq. 
Notice that the appearance of ^Jmq above is only an artifact of the change of variables; 
written in terms of cr, Eq.(20) contains no such term. The theoretical ratios of the first 
two Laplace moments Eqs.(6) and (20) must now be confronted with a corresponding 
ratio involving the experimental data on the J/ip and the T systems. In the case of the 
J/ip one parametrizes the data by a sum of two narrow resonances below DD threshold, 
followed by a hadronic continuum modelled by perturbative QCD. In the case of the T, 
three narrow resonances below the BB threshold are required, at least in principle. One 
obtains 

3 1 1 i'^ 

n(a)Uxp = 2-^^r^'Mvexp(-cTM2) + - / ds exp(-as) Im Ii{s)\QCD ■ (21) 

The experimental ratio is then calculated using Eq.(21) in Eq.(5). The continuum thresh- 
old So is chosen at or below the DD (BB) threshold. Reasonable changes in the value of 
So have essentially no impact on the results, as II{(t) is saturated almost entirely by the 
first two J/V' narrow resonances in the case of charm, and the first T state in the case 
of beauty. In the theoretical ratios the following current values of the QCD parameters 
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have been used: A = 200 - 300 MeV, for four flavours, and A = 100 - 200 MeV, for five 
flavours [|, and < a.G^ >= 0.063 - 0.19 GeV^ 0. 



In the case of charm, one finds that theoretical and experimental ratios match in the 
wide sum rule window: a ~ 0.8 — 1.5 GeV~^, for rUc = 1.39 — 1.46 GeV in the fully 
relativistic case, and a ~ 0.6 - 0.8 GeV^^ = 1.40 - 1.53 GeV in the non-relativistic 
case. For values of a inside the sum rule window, the hierarchy of the various terms in 
the non-relativistic Laplace ratio (20) guarantees a fast convergence. In fact, the leading 
correction in l/m^ is at the 15-20% level, the radiative correction and the non-perturbative 
contribution amount both to less than 10% . At the same time, the next, and next-to- 
next to leading (in l/rric) terms everywhere in Eq.(20) are safely small, as it can be easily 
verified from Eq.(20) noticing that if a 1/2 GeV^^, then r ~ 2mc. Clearly, the complete 
analysis at the level of accuracy of these next-to-leading mass corrections would require 
the evaluation of the perturbative O(a^) terms. Combining the results from both versions 
of the Laplace ratios, leads to the result 

m^(Q2 = = i_4g ± 0.07 GeV . (22) 

In the case of beauty, at small and intermediate values of a the T(IS') provides the 
bulk of the hadronic contribution, i.e. the T(25'), T(35'), and the continuum represent 
a small correction, below the spread in the theoretical ratio due to variations in A and 
in < asG"^ >. Theoretical and experimental ratios match inside the wide regions: cr ~ 
0.4 - 0.8 GeV^^ for m,, = 4.63 - 4.67 GeV in the fully relativistic case, and a ~ 0.20 - 
0.35 GeV~^, for = 4.69 — 4.77 GeV in the non-relativistic one. In the latter case 
all correction terms in Eq.(20) are at the safe level of a few percent. The subleading 
quark mass corrections, though small, are important. For instance, the term of order 
0{as/rnby/rnb) in Eq.(20) is of the same size and sign as the non-perturbative term. 
Hence, it is not fully justified to keep the latter and ignore the former. After combining 
the results from both methods one predicts 

rribiQ'^ = ml) = 4.70 ± 0.07 GeV . (23) 

When comparing the results reported here, Eqs.(22)-(23), with previous determinations 
based on various versions of QCD sum rules 0, - 0, it is important to know which 
values of A and < asG"^ > have been used, as well as which renormalization point has been 
chosen, e.g. some authors determine mQ{Q'^ = —ttlq), which is related to the on-shell 
mass mQ{mQ) through 

mliml) = mli-mDil + ^^"s) • (24) 

The results from the present method are in very good agreement with those of |0| and 
Comparison with other analyses is often made difficult by the lack of information on 
the specific values used for A and the gluon condensate. 
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